Quarantine generated phase transition in epidemic spreading 
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Abstract 

We study the critical effect of quarantine on the propagation of epidemics on an adaptive network 
of social contacts. For this purpose, we analyze the susceptible-infected-recovered (SIR) model in 
the presence of quarantine, where susceptible individuals protect themselves by disconnecting their 
links to infected neighbors with probability w, and reconnecting them to other susceptible indi- 
viduals chosen at random. Starting from a single infected individual, we show by an analytical 
approach and simulations that there is a phase transition at a critical rewiring (quarantine) thresh- 
old w c separating a phase (w < w c ) where the disease reaches a large fraction of the population, 
from a phase (w > w c ) where the disease does not spread out. We find that in our model the 
topology of the network strongly affects the size of the propagation, and that w c increases with 
the mean degree and heterogeneity of the network. We also find that w c is reduced if we perform 
a preferential rewiring, in which the rewiring probability is proportional to the degree of infected 
nodes. 

PACS numbers: 64.60.Aq, 87.10.Mn, 89.75.-k 
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I. INTRODUCTION 



The representation of interactions in a human society as a complex network, where nodes 
and links play the respective roles of individuals and their contacts, has been useful for 
modeling, studying, and understanding many problems in epidemiology l|, |2J. Usually it is 
assumed that diseases evolve faster than the topological evolution of the underlying network, 
so that the links of the network can be regarded as static. With modern mass media, however, 
the presence of an epidemic can be broadcast much faster than disease propagation. This 
information will inevitably change the behavior of the individuals comprising this network 
as, for example, they try to avoid contacts with infected people. In this way a feedback loop 
between the state of individuals and the topology of the network is formed. Networks that 



exhibit such feedback are called adaptive or coevolutionary networks 



Public health services are constantly searching for new ways to try to reduce the spread of 
diseases. Interventions like vaccination js] or total closure of workplaces and schools are very 
effective, but come with a high economic cost. As a less expensive alternative we examine 
here the effectiveness of a "quarantine" strategy, where healthy people are "advised" to avoid 
contacts with individuals that carry the disease. That is, a healthy person has a chance to 
suppress a contact with an infected neighbor and form a new tie with another healthy peer 
(rewiring). The value of this rewiring probability could depend, for instance, on the concern 
that the society has about the disease and, as suggested above, in a globalized world this 
concern will depend on the broadcast news. The degree of media attention about the disease 
is a parameter that could be controlled by public health services. Indeed, it is known that 
spontaneous quarantine in the recent H1N1 pandemic was found to have a large impact in 
reducing the final size of the epidemic 

Based on these observations, we propose two strategies for the propagation of epidemics 
on an adaptive network, in which individuals alter their local neighborhoods with constant 
quarantine probability w as described above, and systematically study the effect that the 
quarantine has on epidemic spreading. We find a phase transition at a critical threshold 
w c above which the epidemic is stopped from spreading. We show also how the epidemic 
spreading in the presence of quarantine depends on the contagion and recovery parameters. 
More importantly, we find that the initial structure of the network plays an essential role in 



disease propagation at criticality, unlike in previous related models 



121 ] where results 



only depend on the average network connectivity. We also introduce a generalized form of 
quarantine, where the probability of rewiring is proportional to the degree of the infected 
nodes, producing a more efficient isolation of the nodes with high degree. This preferential 
rewiring is more efficient than the case with w constant. 



II. ANALYTICAL APPROACH 



We consider the susceptible-infected-recovered (SIR) epidemic spreading model, which is 
well established and accurately describes diseases such as seasonal influenza, SARS, or AIDS 
0, 13)]. Initially, all nodes are in the susceptible state (s), with one node chosen at random 
(seed) in the infected state (i). In the first strategy, strategy A, at each time unit, every 
infected node i in the network tries to transmit the disease to each susceptible neighbor s 
with infection probability 0. If s does not get infected, then with rewiring probability w it 
disconnects its link from % and reconnects it to another randomly chosen susceptible node, 
different from its present neighbors. Thus, the rewiring probability w measures how fast 
susceptible nodes react to the disease (the quarantine probability). Infected nodes recover 
(r) after a fixed recovery time tn since they first became infected, remaining in the recovered 
state forever. 

To estimate w c we start by assuming that the network has a tree structure [14J], so that 
the disease spreads out from the seed and reaches a susceptible node s through only one of its 
neighbors i. Then, if i becomes infected at time to, the probability that s becomes infected 
by % at time t + n, with n = 1,2, .Jr, is (3(1 — /3) n ~ 1 (l — w) 71 " 1 . This is the probability 
that s has neither become infected nor disconnected from i up to time to + n — 1, times the 
probability that i succeeds in transmitting the disease to s at time to + n. Therefore, the 
overall probability that s becomes infected before % recovers is given by the sum 



71=1 

0(1 -[(1-0X1 _«,)]** 



l-(l-/3)(l-w) 



(1) 



This expression for the transmissibility T^ w is equivalent to the corresponding expression 
in the standard SIR model Y^n=\ R — fi) n ~ l [l], but with a non-infection probability 
(1 — — w), instead of (1 — (3). When w = 0, Eq. ([T|) reproduces the known value for 



the transmissibility T = 1 — (1 — f3) tR for the SIR model on static networks In this 
formulation w plays the role of a control parameter of the transmissibility: for fixed values 
of tR and (3, T^ w can be reduced by increasing w. By reducing T^ w we can go from a regime 
in which the epidemic spreads over the population (epidemic phase) to another regime where 
the disease cannot spread (disease-free phase). The transition from the disease-free phase to 
the epidemic phase corresponds to the average number of secondary infections per infected 
node becoming larger than one, allowing the long-term survival of the disease and thus 
ensuring the epidemic spreads to a large fraction of the population. In our problem, the 
expected number of susceptible neighbors that a node has when it just becomes infected is 
given by k — 1, where k — 1 = (k 2 )/{k) — 1 is called the branching factor, and (k) and (k 2 ) 
are the first and second moments, respectively, of the degree distribution P(k). Since T^ w 
is the overall probability to infect a neighbor, the mean number of secondary infections per 
infected node is 



Nf(w) 



(2) 



The infection will die out if each infected node does not spawn on average at least one 
replacement, so for a very large system, the critical point is given by the relation (k — 



1, or 



^-i)p{i-[{i-m-^)] tR } _ 



(3) 



The transition between free-disease phase and epidemic phase is analogous to the static 



link percolation problem 15] , in which each link in a network is occupied with probability 
p and empty with probability q — (1 — p). When p becomes smaller than a percolation 
threshold p c , the giant connected component disappears. In general, p c depends on the size 



of the network N 



161 ]. but in the thermodynamic limit it can be expressed as p c — 1/(k — 1) 



17(. Identifying p with the transmissibility Tg , and using the relation between p c and k 



we find that on the epidemic/disease- free transition line: 



rpA 



Pc- 



(4) 



This result shows that the transition point depends on the initial topology of the network 
through the moments of the degree distribution, as we shall confirm via simulation. 
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A better strategy would be try to avoid the contact between susceptible and infected 
individuals before the attempt to infection, i.e., only avoid the contact when you know an 
individual is sick. In this second strategy, strategy B, at each time unit every susceptible node 
attached to an infected node disconnects its link from i with probability w and reconnects it 
to another randomly chosen susceptible node, different from its present neighbors. If s does 
not rewire its link, the infected node % tries to infect it with infection probability (3. Notice 
that Nf(w) = (1 — w)Nf(w). For this strategy Eq.flTJ is replaced by 



g 0(1 -«;){l- «;)]** ,} 
T ^ w = l-(l-/3)(l-u;) • (5) 

Comparing Eq.(JTJ and Eq.fjS]) we see that T^ w < T^ w for identical values of (3 and w. The 
new condition for disease to die in strategy B out follows from Eq.(j5]) 

(« _ 1)0(1 -w c )h- [(i - 0)(i - w c )] tR \ 

= 1. (6) 

1-(1-0)(1-™ C ) {) 

The two strategies represent different scenarios depending on the knowledge of the state 
of infection of the nodes. For strategy A, susceptible nodes have no information about 
the state of their neighbors until they are in physical contact. In contrast, for strategy B 
susceptible nodes know the state of their neighbors before they are in physical contact. We 
will show that this difference in the knowledge of the states of the nodes results in strategy 
B being more effective at stopping epidemic spreading. 

Fig. [T] shows phase diagrams in the w — (3 plane for strategy A and strategy B, obtained 
by the numerical solutions of Eq. ((3]) and Eq. (]6]) respectively, for a network with k = 5 
(p c = 0.25). Two phases emerge: the epidemic phase for T^f > p c , and the disease-free 
phase for T^f < p c . The location of the critical line separating the two phases depends on 
the recovery time £r. Fig.QJa) illustrates the two limiting cases for strategy A: tn — 1 (upper 
curve) and = oo (lower curve), which are given by the expressions j3{t^ = 1) = — 1) 
and /3(£r = oo) = w/(n + w — 2). Therefore, the pairs of (w,(3) values in the region above 
the curve (3(tn = 1) are always in the epidemic phase, and below the curve (3(tR = oo) 
are always in the disease-free phase. A striking consequence is that if (3 is larger than the 
percolation threshold 1/(k — 1), the propagation of an epidemic cannot be stopped, even 
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FIG. 1. j3-w phase diagram for p c = 0.25. The curves correspond to t# = 1 (upper) and £r ; = oo 
(lower). The dark gray region will always be an epidemic phase, and the light gray always a 
disease- free phase: (a) strategy A, (b) strategy B. 



with the largest rewiring probability w = 1. Thus, strategy A is not an efficient mechanism 
to control an epidemic when /3 is higher than p c . 

For strategy B, Fig. H(b) illustrates the two limiting cases tn = 1 (upper curve) and 
tn = oo (lower curve), which are given by the expressions /3(tn = 1) = 1/(k — 1)(1 — w) and 
/3(tn = oo) = w/(k — nw + 2w — 2). In contrast to strategy A above, if k does not diverge, 
i.e., the network has finite p c , the epidemic can always be stopped, even for — > 1. Notice 
that the maximum value of w needed to stop a epidemic is w c = (k — 2)/(k — I) = I — p c . 

Our theoretical predictions illustrate a novel feature about the dynamics of adaptive 
networks in SIR models. While previous adaptive SIS and SIRS models have transition 
values that depend only on the average connectivity of the network (k), and are independent 
of the heterogeneity or structural correlations of the initial topology 3] , our SIR model 
predicts dependence on the topological structure through the higher order moments of the 
degree distribution. 



III. SIMULATION RESULTS 

Our analytical approach predicts, through Eq. (jHJ) and Eq. (jBJ), a dependence of w c on 
the initial network. In order to test and explore this dependence, we performed extensive 
numerical simulations of our model starting from different topologies. We first used Erdos- 
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Renyi (ER) networks with Poissonian degree distribution P(k) = (k) k jk\. Even though 
these types of homogeneous networks are common in nature, many real social networks 
are well represented by heterogeneous networks. Thus, we also used finite scale-free (SF) 
networks with degree distribution P(k) = k~ x exp(—k/K) /Li\(e~ 1 ^ K ), where K is the degree 
cutoff. This distribution represents networks with a finite threshold p c , and appears in a 
variety of real- world networks fisl. fiol] . We only consider epidemic propagation on the largest 
connected cluster of the network, the giant component (GC). 

In Fig. [2] we compare both strategies for ER networks. We plot as a function of w, the 
average fraction of infected nodes nj(w) = Nj(w)/Ngc, where Nqc is the size of the giant 
component, on a ER network. This fraction is normalized by the corresponding fraction on 
an identical fixed network, i.e., for the SIR model without quarantine, nj(0). We can see 
that w c for strategy B is lower than w c for strategy A, as expected. This relation will hold 
for any topology with the same k (see Eq.([n]) and Eq.([2])). Since strategy B is more effective, 
from here to the end of the paper we will show only simulation results for strategy B. 




FIG. 2. Plot of ni(w)/ni(0) as a function of w for strategy A (0)> an d strategy B (□), for a ER 
network with (k)cc = 4.07, N = 10 4 and 10 4 realizations. As we predicted previously, strategy B 
is better than strategy A, shown by < w^. 

Fig. [3] shows, ni(w)/ni(0) vs w for different (3. In this figure, we can observe the strong 
effect of quarantine and the critical threshold w c above which ni(w) approaches zero. Thus, 
in the disease-free phase (w > w c ) only a small number of individuals get infected and the 
disease quickly dies out. We observe, as expected, that the values of w c increase with 0. 
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This behavior matches the phase diagram of Fig. [T] (b), where the critical line has a positive 
slope. Scaling the horizontal axis by the values of w c obtained by numerically solving Eq. ([6]) 
collapses the curves, showing a excellent agreement between theory and simulations, as well 
as a scaling behavior of the form rij(w) = 71/(0)/ (w/w c ), as shown in the inset of Fig. [3J 




FIG. 3. Fraction of nodes infected as a function of the applied quarantine parameter w divided by 
the fraction of infected nodes at w = for three different infection probabilities j3 in an ER network. 
(O) P = 0.05, (□) /3 = 0.1 and (0) P = 0.15. Above a threshold value of w, no finite fraction of 
the network becomes infected. (Inset) Data with w rescaled by the appropriate w c calculated from 
Eq. ([6]). The curves collapse very well, showing universal behavior and good agreement with the 
theory. In all simulations, {k)cc = 4.07 ((k) =4 over all nodes), N = 10 and averages are over 
10 4 realizations. 

We also explored the dependence of w c on the connectivity of the network, by computing 
ni(w)/ni(0) vs w for different (k) (see Fig. [4]). We observe that w c increases with (k), due 
to the fact that propagation is facilitated by having more neighbors, as in the original SIR 
dynamics. The inset of Fig. H] shows the collapse of all curves. Again, the good agreement 
between theory and simulation confirms that Eq. ([6]) is a valid expression of the transition 
point for the adaptive SIR model on ER networks. 

Given that the second moment, and therefore k, is large in heterogeneous networks, the 
critical value w c turns out to be larger in heterogeneous networks than in homogeneous 
networks with the same mean degree (k) and size N, as seen by considering Eq. ([6l) in the 
tR ^> 1 limit, where 

Wc ^ i 1 ■ (7) 

/3(«-2) + 
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FIG. 4. Fraction of nodes infected as a function of the applied quarantine parameter w divided by 
the fraction of infected nodes at w = for ER networks with three different (k): {k)cc = 4.07 (0)> 
{k)cc = 6.015 (□) , {k)cc = 10 (0)- Again, while below w c a finite fraction of the network become 
infected above w c the epidemic vanishes. (Inset) Data rescaled by w c calculated from Eq([U]) . In 
all simulations, /3 = 0.05, N = 10 4 and averages are over 10 4 realizations. 



Since k for heterogeneous networks is much bigger than in homogeneous networks we 
expect that w c increases as the heterogeneity increases with larger K. For k — > oo we 
expect that w c — > 1 and that the transition will eventually disappear for large heterogeneous 
networks. 



In Fig. |5l we show simulations on scale- free networks for different values of K 20]. In 
good agreement with our predictions from Eq. (JTJ), as k increases with K, w c increases. In 
the inset we rescale by w c obtained from Eq. flfJJ), and find good collapse (Fig. H^inset)), 
confirming the general validity of Eq. (E]) for heterogeneous SF networks. 

In very heterogeneous networks it is well known that due to high degree nodes, propaga- 
tion processes are very difficult to stop 15|, |2l| . Thus, a constant probability of rewiring is 



not effective for controlling epidemics in these networks, because such a strategy assumes all 
nodes have the same importance in an epidemic, ignoring the function of higher degree nodes 
as superspreaders. It is also well known that removing the high degree nodes of the network 



(targeted percolation) 



15 



22 



23| is a more efficient method to stop propagation processes 



than random removal. This type of strategy is expected to be superior, but requires global 
information about the network. 



To test this prediction, we propose a new strategy of type B where w depends on the 




FIG. 5. Fraction of infected nodes still infected as a function of w for (3 = 0.05, tn = 20 on an SF 
network with K = 10 (O) and K = 20 (□), A = 2.1 with N = 10 5 . p c ^ 0.375 for K = 10 and 
p c = 0.215 for if = 20. (Inset) Rescaled by w c . The result shows an excellent agreement between 
the theory and the simulation. 



degree k of the infected node, with Wk given by the general form 



w k = w k (a) = ^k a 



(8) 



where 7 is a constant that controls the highest possible value of w^, and it is equal or smaller 
than kmtxj where k max is the largest degree of the network. For a = and 7 6 [0,1] we 
recover the results for a constant value of w, with 1/7 = 7. For a > and 7 = k^ x the 
rewiring increases with the degree k of the infected node, and decreases with increasing a. 
In the limit of a — > 00 the rewiring process is equivalent to a targeted rewiring where only 
the links of the highest connected node(s) are rewired. To compare the cases with a = 
and a > 0, we use the average (w) over the network, choosing a for the targeted case such 
that (w) is equal to w in the uniform case. 

In Fig. we plot nj(w)/ni(0) as a function of (w) for SF networks with K = 10 and 
= 20. As expected, w c for a > is lower than for a = 0. The preferential rewiring 
reduces the value of w c because quarantine is more effective at isolating the superspreaders. 



IV. SUMMARY AND CONCLUSIONS 



We have introduced and studied two strategies for the propagation of epidemics on evolv- 
ing networks of social contacts. The states of the individuals are changed according to SIR 
dynamics, while the network evolves according to a quarantine mechanism based on lo- 
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FIG. 6. Fraction of infected nodes still infected as a function of (w) for /3 = 0.05, tR = 20 on an 
SF network, A = 2.1 , N = 10 4 and K = 10, 20. For general strategy B with a = 0, (K = 10 (■) 
and K = 20 (A)), and with a > 0, (X = 10 (□) and .ff = 20 (a), each (w) represents a different 
a. 



cal information, in which susceptible individuals replace their infected neighbors by other 
susceptible peers with probability w. We demonstrated by an analytical approach, and con- 
firmed by numerical simulations, that the size of the epidemics can be largely reduced by 
increasing the probability of rewiring, and that the propagation can be eventually stopped 
by using high enough values of w. In other words, quarantine is an effective way to halt the 
appearance of epidemics that would otherwise emerge in the case of a static network. For 
strategy A, when the infection probabliy is larger than a threshold the quarantine mechanism 
is not effective any more and disease propagation becomes unavoidable. This is because the 
quarantine model only breaks contact after allowing for a chance of infection, thus for high 
enough infection probability, the spreading is irreversible, even with a rewiring probability 
of unity. For strategy B, the quarantine mechanism is effective in homogeneous networks 
and finite heterogeneous networks, even for infection probability equal to unity. The tran- 
sition disappears only for large, very heterogeneous networks. In these scale-free networks, 
the critical rewiring probability becomes very large and the transition eventually disappears 
for large enough systems, thus only the epidemic phase is observed. This is likely due to 
the presence of individuals with very large connectivity that can spread the disease over a 
large fraction of the population, even for small infection probabilities. To confirm this, we 
introduced a generalized form of strategy B, where the quarantine depends on the degree 
of the infected nodes. This preferential rewiring isolates the superspreaders more efficiently, 
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reducing w c and preserving a finite transition even for a scale-free networks. Lastly, in SIR 
dynamics the final frozen state, where everybody is either recovered or susceptible, is reached 
rather quickly — in a time of order IniV according to our simulations — and thus an SIR 
network does not evolve much, so the initial topology is preserved during the entire evolu- 
tion of the system. In contrast, in adaptive models with SIS or SIRS dynamics the system 
evolves for a very long time in the epidemic phase — times that grow exponentially larger 
with iV — thus the network moves towards a stationary topology similar to an Erdds-Renyi 
network, independent of the initial topology, and rendering the initial topology irrelevant. 
Therefore, unlike other adaptive network models of epidemic spreading, in our model the 
epidemic threshold has a strong correlation with the topology of the network, which remains 
relatively unchanged at criticality, and the social structure of connections when epidemics 
begin to propagate is crucial in the state of the final outbreak. 
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